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ABSTRACT 

We study the prospects of finding the first quasars in the universe with ALMA and JWST. 
For this purpose, we derive a model for the high-redshift black hole population based on observed 
relations between the black hole mass and the host galaxy. We re-address previous constraints 
from the X-ray background with particular focus on black hole luminosities below the Eddington 
limit as observed in many local AGN. For such luminosities, up to 20% of high-redshift black 
holes can be active quasars. We then discuss the observables of high-redshift black holes for 
ALMA and JWST by adopting NGC 1068 as a reference system. We calculate the expected 
flux of different fine-structure lines for a similar system at higher redshift, and provide further 
predictions for high-J CO lines. We discuss the expected fluxes from stellar light, the AGN 
continuum and the Lyman a line for JWST. Line fluxes observed with ALMA can be used to 
derive detailed properties of high-redshift sources. We suggest two observational strategies to find 
potential AGN at high redshift and estimate the expected number of sources, which is between 
1 — 10 for ALMA with a field of view of ~ (1')^ searching for line emission and 100 — 1000 for 
JWST with a field of view of (2.16')^ searching for continuum radiation. We find that both 
telescopes can probe high-redshift quasars down to redshift 10 and beyond, and therefore truely 
detect the first quasars in the universe. 
Introduction 



After the firs t detection of su permassive black 
holes at z ~ 6 ( Fan et al.l 12001 ). there has been 



a strong effort to find more sources at high 
redshift and to study their properties. Fur- 
ther high-z quasi-stellar objects (QSOs) have 



Survey fSDSS) 


(Fan et al.l 2001; Carballo et al. 


2006': 'Cool et al. 


20061 


Fan et all 2006; Shen et al. 


2007; Jnada et al. 2008 


), the Canada- France high- 



trosco pv (I Juarez et al. 20081). mid-IR observa- 
tions ( Stern et al. 2007 ; iMartmez-Sansigre et al.l 
2008h . Swift observatio ns tSambruna ct af. 20Q7p , 
VLB A observations (jMomjian fc Carillil 120071 ). 
the Very Large Array as well as the Very Long 
Baseline Array (Momii an et al. 20081 ). Some of 
their host galaxies have als o been detec t ed in 
CO and submm emission ("Omont et al.' 1 19961 



Carilh et al.1 ^2002: Walter et al. 2004; WeiB et al 



2007; Riechers et al. 2008 



UKIDSS Large Area Survey (IVenemansI 120071) 



and the Palomar-QUEST survey (jBauer et al 



20061 ). They have been st udied by follow-up ob - 
servations with Chandra ("Shemm er et al. l2006h . 
the Subaru Telescope (,GotQ ,2007t) . with IR 
spectroscopy ( Juarez et al.l 20071) . near-IR spec- 



These observations furthermore inferred a num- 
ber of relations between the mass of the cen- 
tral black hole and the properties of their host 
galaxy. The Magorrian relation connects the 
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bulge with the black hole mass ( 


Magorrian et al. 


1998 


; iFerrarese fc Merritt' '2000' 'Cebhardt et al. 


2000 


; Graham et al. 


200L .Merritt & Ferrarese 


2001 


; Tremaine et al.l 


2002; Haring & Rix 2004 


Petersonl 20081 Somerville et al.l 20081). Quasars 



with black hole masses of the order 10^ Mq are 
observ ed with supersolar m etallicities even at 
z ~ 6 ([Pentericci et al.l 120021 ). which is indicated 



1 



in further studies as well (iFreudling et al. 2003t 



Maiolino et al. 2003: Dietrich et al 



2008al lcl) 



In 



the nearby universe, correlations have b een found 
betw e en black hole mass a nd metallicity (jWarner et al 
20031: iKisaka et all l2008l ). with a mild slope of 
0.38 ± 0.07. This is plausible, as it is gen- 
erally recognized that the overall mass of a 



galax y is correlated with i ts metallicitv (IFaber 
19731: [Zar itskv ct al. 199^ Ijablonka et all Il996l : 



Trager et al.ll2000f) . In this paper, we will in gen- 



eral assume that the black hole mass is a good indi- 
cator for the evolutionary stage of the quasar host 
galaxy, in particular with respect to the Magorrian 
relation and the expected metallicity. 

In the local universe, active galactic nuclei 
(AGN) can be studied in much more detail 
and with higher resolution. An example that 
has been studied particularly well, in differ- 
ent wavelengths and with high resolution is the 
AGN NGC 1068. Its direct X-ray emission 
is shielded through absorption along the line 
of sight, but there is evidence for reflected X- 
ray p hotons that are scattered into the line of 

Its molecular 



sight (jPounds fc Vaughan 200' 



disk has been studied through emission from 
different species, in particular molecular CO 



and fine-structure lines ( 


Schinnerer et al.l 


200C 


GalUano et al. 2003, 2005 




Soinoglio et al.l 


2005 


Davies et all l2007bllal: IPoncelet et all 2007f). 


The 



continuum emission and the distribu tion of dust 
has b e en measured i n furt h er studies (iKrips et al 
20061: 



Mason et al.'' '2006": 



'Rhee & Larkin '200? 



Tomono e t al. 2006; Howell ct al. 2007; Honig ct al 
20081 ). Such observations in particular allow 
to probe stellar and gaseous dyn amics and the 



galac ti c gravitational poten ti al (iGerssen et al 



20061: lEmseUem et"aLl (20061: iDas et al.l l2007l) 



Observations with the Green Bank Telescope 
(GBT) provide evid e nce for water maser activit:v 



((jti 1 ) provide evia e nce tor water maser activity 
( Greenhill fc GwinnI Il995l : iGreenhill etHI Il996l : 
2006f l. The inner radio jet can 



Kondratko et al 



be detected with ob servations at 7 mm wavelength 
(jCotton et al.ll2008l ) . Given this amount of detail, 
NGC 1068 thus provides an excellent test case for 
the inner structure of active galaxies. 

The formation of quasars and their supermas- 
sive black holes is still one of the unresolved rid- 
dles of structure formation and cosmology. The 
simplest scenarios assume that they have grown 
from the remnants of the first stars, which are 



believed to be very massive jAbel et al. 2002 : 
Bromm fc Larsoiil 2004 ; Glover 2005 ) , and whose 
black hole remnants could grow further by accre- 
tion. Such a scenario however has problems, as 
the remnants of the first stars typically do not end 
up in the most massive quasars at redshift z ^ 6 
( Trenti fc Stiavellil 2008), and radiative feedback 
from the stellar progenitor ca n delay accretfon as 
well ( Johnson fc Bromm 2007: 1 Alvarez et al.ll2008t 
Milosavlievic et al.ll2008l ). In case of Eddington 
accretion, seed black holes of ~ 10^ Mq are re- 
quired in order to grow to the observed supermas- 
sive black holes at z ~ 6 (IShaDiroll2005h . 

Recently, it has also been discussed whether 
the first stars in the early universe were powered 
by dar k matter annihilation rather than nuclear 
fusion (jSpolvar et all liooj lloccd '2008V Such 
stars could reach masses of the order 1000 Mq 



(jPreese et al.l l2008t llocco et al.l l2008l ) and were 



considered as possible progenitors for the first su- 
permassive black holes. The evolution of such 
stars on the r tiain s eque nce has been calcul ated by 
Taoso et al.l (|2008h and lYoon et"aD (|2008h . How- 
ever, it was shown that such stellar models are 
highly constrai ned by the observed re ionization 



optical depth (jSchleicher et al.l l2008bl lah. Also, 



we note that such seeds would still require super- 
Eddington accretion to grow to the observed su- 
permassive black holes at z ~ 6. 

Therefore, alternative scenarios have been con- 
sidered that lead to the format fon of massive seed 
black holes by direct collapse jEisenstein fc Loeb 
1995; Koushiappas et al. 2004 : iBegelman et al 



2006[ ISpaans fc Silkl 120061 : IPiikstra et all l2008h . 

However, for such scenarios, the situation is also 



not en tirely clear. For instance. iLodato fc Nataraian 
( 20061 ) argue that the gas in such hal os should frag 



ment if H2 cooling is effici e nt, an d lOmukai et al 



(|2008l ) and iJappsen et al.l (|2008l ) suggest that a 
non-zero metallicity will lead to fragmen tation as 
well. Simulations by IClark et al.l (j2008r ) support 
these conjectures. 

The black hole population at high redshift is 
further constrained by observations of the soft X- 
ray background. These constraints indicate that 
the population of high-redshift q uasars was not 



suffic ient to reionize the universe (jPiikstra et al 



20041 ). and suggest an upper limit to the black hole 
densi ty of 4 x 10" Mq Mpc"^ (jSalvaterra et al. 
20051 ). The latter depends in particular on the 



2 



adopted spectra and luminosity of the black hole, 
as well as their duty cycles, and it has been argued 
that lar ger black hole p opulations are conceivable 



tnat lar ger black nolo p opulations are conceivable 
as well ( Zaroubi et al.l I It is therefore im- 

portant to probe the number density of black holes 
and AGN activity as a function of redshift by di- 
rect observations. 

In this work, we explore the possibility of 
searching for the first quasars in the universe 
with ALMA and JWST B ALMA is par- 
ticularly relevant in light of recent predictions 
on the observability of the first miniquasars 
(jSpaans fc MeiierinkI 120081) and high-re dshift Ly- 
man a galaxies (jPinkelstein et aL |2008[ ) , while it 
is one of the main scientific goals of JWST to 
find the first galaxies at high redshift. In [J2l we 
introduce a model for the black hole population 
which is based on observed properties and corre- 
lations between black holes and their host galaxy. 
In 331 we discuss how the number of active black 
holes is constrained by the unresolved X-ray back- 
ground. In 21 we analyze the observed properties 
of NGC 1068 and derive further predictions for 
very high- J CO lines. In !j5l we use these in- 
formation to derive the relevant observables for 
ALMA and JWST and predict the number of de- 
tectable sources as a function of redshift. Further 
discussion and outlook is provided in 

2. The high-redshift black hole population 

We model the population of high-redshift black 
holes in the framewor k of the Sheth-Tormen 

'2001'), a modified ver- 

Press fc Schechteil (|l974[ ) method. 



formalism (|Sheth et a l. 
sion of the 

Both prescriptions agree reasonably well with 
simulations of high-redshift structure formation 
(jjang-Condell fc Hernguist 2001 : Heitmann et al 



-nqi 

20061 : iReed et al.l 120071: lOreif et al.l l2008l) . This 



formalism predicts the comoving number density 
of halos with different masses, dn/ dXogMh, where 
Mh is the halo mass. We calculate the virial radius 
Tvir assuming a mean overdensity of 178 and esti- 
mate the velocity dispersion a as the Kepler veloc- 
ity at this radius. For simple profiles such as the 
singular isothermal sphere, the velocity dispersion 
is independent of position, while in more realistic 
scenarios, there will be some scatter around this 



^ http: / /www. eso.org/ sci / facilities / alma/index. html 
^http: / /www. stsci.edu/jwst/overview/ 
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Fig. 1. — The co-moving number density of halos 
that contain black holes as a function of black hole 
mass for different redshifts. 



value. We adopt the correlation between black 
hole mass Mbh and ve locity dispersion a given by 
(|2002l) as 



Tremaine et al 



log 



BH 



/31og( - 



(1) 



with a = 8.13 ± 0.06, /3 = 4.02 ± 0.32 and cro = 
200 km/s. The expected black hole population is 
then given in Fig. [1] This shows that a significant 
number of potential sources is available down to 
z ~ 12, which however drops significantly towards 
higher redshifts. 

If black holes accrete efficiently and release the 
gravitational energy as radiation, the luminosity 
can be estimated with the Eddington luminosity 
as 

'Mbh^ 



L 



Edd 



= 3.3 X 10'* 



M, 







Lr. 



(2) 



On the other hand, accretion may also be lim- 
ited by gas dynamics, and the Bondi solution for 
spherical symmetric accretion shows that the ac- 
tual accretion rate may be considerably smaller if 
the gas is heated by radiative fe edback or super- 
nova e xplosions. As discussed bv ljoUev fc Kuncic 
([20081), part of the released gravitational energy 
may go into mechanical feedback like jets and out- 
flows rather than radiation. This is consistent 
with measurements of the black hole mass - lu- 
minosity relation bv feaspi et al. ( 2000l ). which in- 
dicate that typical luminosities are in fact sub- 
Eddington. They derive two correlations between 
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black hole mass and luminosity, depending on 
whether the velocity estimate is obtained taking 
the mean (case A) or the root mean square (rms, 
case B). They use the parametrization 

Mbh = a X 10^ X ^ Mq. 3 

\ 10^^ erg s ^ / 

In case A, they obtain qa = ^^'^ = 

0.545 ± 0.036, in case B, they find an = 5 .75t°il 
and bs = 0.402 ± 0.034. Like lKaspi erall lIpOO), 
we adopt a factor of 9 to correct for the bolo- 
metric luminosity. In both cases, the luminosity 
scales roughly with and most quasars are 

sub-Eddington. As shown in the next section, 
the X-ray background constraints on the black 
hole population are much less significant with this 
mass-luminosity relation. 

3. Constraining the relative number of ac- 
tive black holes 

The population of high-redshift black holes 
is constrained by the unresolved soft X-ray 
background (SXRB). The total SXRB in the 
energy range 0.5 — 2 keV is (7.53 ± 0.35) x 
10~^^ erg cm~^ s~^ deg~^, with 94 ^? per cent 



made up of discrete X-ray sources fjMoretti et al 
20031 ). The unresolved fraction of the SXRB 



was re-analyzed by iDijkstra et al.l (|2004f ). pro- 



viding a mean and maximum intensity of the un- 
resolved SXRB flux of 0.35 x lO^^^ g^^^ ^ 23 x 
10"-^^ erg cm~^ s~^ deg'^, respectively. This 
corresponds to specific fluxes per energy interval 
of the order 0.7 - 2.5 keV cm^^ g-i gj.-i keV-^ 
at 1 keV, which is the upper limit for possible 
contributions from the early universe. 

3.1. Assumptions 

For the X-ray sources, we coii s ider s pectra 
of similar type as ISalvaterra et al. I (|2005l ). con- 
sisting of a soft multicolor disk (MCD) compo- 
nent and a power-law ( PL) c omponent at high 
energies (iMitsuda et all Il984l: iMiller et al.l l2003t 
Miller fc ColbertI l2004[ ). The MCD component 
peaks at frequencies 



hv - 0.095 keV 



BH 



106 Mq 



-1/4 



(4) 
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Fig. 2. — The expected contribution to the SXRB 
in keV cm~^ sr~^ keV~^ as a function of the 
active quasar fraction fa and the spectral param- 
eter $, under different assumptions for the typ- 
ical black hole luminosity. For a spectrum with 
the same amount of flux in the MCD and the PL 
component, we expect $ < 0.5. Upper left: Ed- 
dington luminosity, black holes in the mass range 
from 1 to 10^ Mq. Upper right: Eddington lumi- 
nosity, black hole masses in the range from 10^ to 
10^ Mq. Lo wer left: Case A b lack hole - luminos- 
ity relation ( Kaspi et al.ll2000l ). Lowe r right: Case 
B bl ack hole - luminosity relation (jKaspi et al 
20001 ). For the latter models, results are indepen- 



dent of the lower mass cutoff due to the scaling of 
the employed luminosity relations with mass. In 
general, we find that models assuming Eddington 
luminosity find significantly stronger constraints 
than models following the observed luminosity re- 
lation. This implies that the fraction of active 
quasars may be significantly higher than for Ed- 
dington luminosities. 



and the spectral energy distribution (SED) scales 
with j/^/"^ for smaller frequencies and rolls off ex- 
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Fig. 3. — An illustration for the inner structure of 
NGC 1068. It shows the central black hole with 
the inner accretion disk radiating in X-rays, as well 
as the inner torus embedded in a large molecular 
disk consisting of individual high-density clouds. 
The jet is sketched as well. The torus is expexted 
to have a size of ~ 3 pc, while the molecular disk 
extends out to ~ 100 pc. In the XDR model con- 
sidered here, we focus on those clouds that are un- 
shielded by the torus and where emission is stim- 
ulated by the impinging X-ray flux. 

ponentially for larger frequencies. The MCD com- 
ponent of high-redshift quasars is therefore red- 
shifted to even smaller frequencies and cannot con- 
tribute to the unresolved SXRB. For the PL SED, 
we assume a scaling according to with a ~ 1. 
It is the PL component of these sources which is 
redshifted into the SXRB. The spectral range in 
which these sources emit is not completely clear. 
Here, we adopt a spectral range from 1 eV to 
500 keV. The ratio r of the flux in the MCD com- 
ponent to the flux in the PL co mponent is ex 



pected to be roughly o f order 1 (jMitsuda et al 



19841: iMiller et al.ll2003l ). but we treat it as a free 
parameter in this analysis. A further fraction fg 
of the radiation may be shielded by the inner 
torus. This fraction is geometry dependent, but 
may be of the order of 20 — 30%. The amount 
of luminosity that can potentially contribute to 
the SXRB thus scales with the effective parame- 
ter $ = (1 — /s)/(l + r). For typical spectra where 
the amount of flux in the MCD and the PL compo- 
nent is comparable, we therefore expect $ < 0.5. 
We further introduce the parameter fa, which de- 
scribes the fraction of active quasars at a given 
redshift. 



3.2. Formalism 

The contribution to the background intensity 
is then given from an integration along the line of 
sight as 

'"'An J H(z){l + zY ' 

where v is the observed frequency, H{z) the ex- 
pansion rate at redshift z and Pu{{l + z)^, z) the 
mean proper volume emissivity for photons with 
frequency v{l + z) at redshift z (in units of energy 
per unit time, unit frequency and unit proper vol- 
ume). To evaluate that, we first calculate the num- 
ber of photons No of frequency v emitted at red- 
shift z within 1 Mpc'^ (comoving), which is given 
as 



M„ 



dn 



dlogMi 



-/aA^(i>,MBH)dl0gM, 



(6) 

where Mmin and M^ax are the (redshift-dependent) 
minimum and maximum halo masses that yield 
black hole masses in a reasonable range which we 
specify below in more detail, and N{i>, Mbh) is 
the number of photons of frequency v emitted per 
second by a black hole of mass A/bh- The number 
of baryons Ni, in 1 Mpc'^ (comoving) is 



Nh = 



Pb 1 Mpc'' 
1.2m„ 



(7) 



where pt is the comoving baryon density and rrip 
the proton mass. The emissivity is then 



P. 



hi>Ntot,i> Pb(l + 



1.2mr, 



(8) 



3.3. Constraints 



With the method described above, we calculate 
constraints for the fraction of active quasars fa as 
well as the spectral parameter $ for four different 
scenarios (see Fig. [2]). The first scenario assumes 
that black holes populate the mass range from 1 
to 10^ Mq and radiate at Eddington luminosity. 
This may be the case if a sufficient number of black 
holes is provided as remnants from the first stars. 
In this case, one finds rather tight constraints that 
do not allow high values for /„ and $ simultane- 
ously. 
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However, as the initial mass function (IMF) in 
the early universe is not yet understood, it is un- 
clear whether such stellar mass black holes are 
present in large numbers. As discussed in the in- 
troduction, supermassive black holes may in fact 
have more massive progenitors that form through 
direct gas collapse in larger halos. We assume here 
a lower cutoff for the black hole mass of ~ 10^ M©. 
That allows for an increase in the fraction of active 
quasars by a factor of a few. 

Finally, as discussed in Sj2l the Eddington lumi- 
nosity is only an upper limit to the actual black 
hole luminosity. It depends on the gas supply as 
well as the amount of energy going into mechan- 
ical feedback. To address alternative scenarios as 
well, we calculate the constraints assuming the 
black hole mass - luminosity relations given by 



Kaspi et al 



( 2000f ). As a consequence, a higher 
fraction of active quasars is consistent with the 
observed background for parameters $ between 
0.1 — 1. Due to the steepness of the luminos- 
ity function, the results are independent of the 
adopted lower cutoff for the black hole mass. In 
case A, fa ~ 50% seems feasible for $ ^ 0.2, while 
in case B, fa ^ 50% requires $ ~ 0.1. For a spec- 
trum with equal amount of flux in the MCD and 
the PL component ((f> < 0.5), active quasar frac- 
tions of 10 — 30% can be adopted. 

With respect to the observational search for the 
first quasars at high redshift, the two latter scenar- 
ios are particularly favorable, as they are consis- 
tent with a considerably larger number of active 
quasars in the early universe. The main caveat 
here is the question whether the observed relation 
between black hole mass and luminosity can in- 
deed be used at these early times. So far quasars 
have always appeared as highly evolved systems 
at any redshift observed, which may imply that 
this relation indeed holds. In any case, the first 
scenario assuming Eddington luminosity provides 
a firm lower limit of about 1% for the fraction of 
active quasars if $ is assumed to be of order 1. 

4. NGC 1068: A case study 

In this section, we review known observational 
properties of NGC 1068, derive predictions that 
will be relevant for high-redshift observations with 
ALMA and JWST and generalize them to objects 
with different X-ray fluxes. 
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Fig. 4. — A model for the X-ray chemistry in 
NGC 1068. The adopted flux impinging on the 
cloud is 170 erg cm~^. The adopted density 
is 10^ cm~^. Top: The abundances of different 
species as a function of column density. Middle: 
The low- J CO lines as a function of column den- 
sity. Bottom: The high-J CO lines as a function 
of column density. 
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4.1. Observed properties 

As discussed in the Introduction, the nearby 
AGN NGC 1068, a Seyfert 2 galaxy at a distance 
of 14.4 Mpc, has been studied in great detail and 
with high resolution. As the wealth of available 
data could easily give rise to a review of its own, we 
restrict the discussion here to only those properties 
that will be relevant with respect to observations 
at high-redshift with ALMA and JWST. 

The mass of the central black hole has been 
estimated from VLBI water maser observations 
to be ~ 10*" Mp) with a luminosity of 50% 



of the Eddington hmit I Greenhill et al. 1996h . 



Even though the X-ray emission from the cen- 
tral black hole is shielded by an absorber on the 
line of sight. X-rays that are scattered into the 
line of sight have been d etected unambiguously 
( Pounds fc Vaughan 2006h . Observations of rota- 
tional CO lines and the H2 2.12 fim line indi- 
cated an inner structure consisting of an extended 
molecular disk and a toroidal structure which com- 
pletely absorbs the X-ray emission along the line 
of si ght. It is slightly til ted with respect to the 
disk fGalliano et al. '2003"). The H2 emission map 
of [Galli ano fc A Uoin (2002) shows a bright knot 
70 pc in the east of the central engine, with a 
size of ~ 35 pc. Th e knot is also visible in CO 
(iGalliano et al ] l2003h . As they show, the observed 
fluxes ca n be explained with t he X-ray excitation 
model of lMalonev et al. I (|l996h when the following 
parametrization is adopted: 

• The central engine is a power-law X-ray 
source with spectral slope a = —0.7 and lu- 
minosity of 10'*'* erg s~* in the 1 — 100 keV 
range. 

• We observe emission from molecular clouds 
with a density of 10^ cm~^, a column of 
10^^ cm"'^ at a distance of 70 pc and solar 
metallicity. 

This region of bright emission consists therefore 
of a number of small clouds with densities of 
10^ cm~^, which in total cover a significant part 
of the telescope beam. In addition, there may be 
further contributions from gas at lower densities, 
which emits less efficiently, but fills a larger vol- 
ume. For the X-ray luminosity of NGC 1068, one 
cannot expect CO line emission much closer to the 



black hole, as the higher X-ray flux then would de- 
stroy CO very effectively, unless it is shielded by 
some absorber. The black hole luminosity thus 
regulates the distance to the central source at 
which the CO emission peaks. This situation is 
sketched in Fig. [3] 

Further studies also measured the fluxes from 
other species, including the relevant fine-structure 
lines. We refer here in particular to the data of 
Spinoglio et ail (|2005l ). which are used m ij5.1l to 
calculate the expected fluxes from higher redshift. 
The situation is similar for the observed amount 
of continuum emission. Here we use the data of 



Brauher et all (|2008f ). 

As JWST will see the stellar light from high- 
redshift galaxies, the star formation rate is a 
highly relevant quantity for estimates regard- 
ing the expected flux. Indeed, observations of 
NGC 1068 find a starburst component from a 
circumnuclear ring of ^ 3 kpc in size, with a 
stellar mass of ~ 10® Mq and an age of 5 Myr 



(jSpinogho et al.M2005f ). On scales of a few hun- 
dred parsecs, one finds star formation rates of a 
few times 10 Mq yr~* kpc~^, which can reach 
~ 10 M(?) yr~* kpc~^ i n the inner few ten par- 



secs ([Davies et al.ll2007bl ). This is close to the star 
formation rat e in Eddington-limited s tarbursts as 
suggested by [Thompson et al. ( 2005 ). The to- 
tal star formation rate of NGC 1068 therefore is 
^ 10 Mq yr-1. 

4.2. Predicted properties 

High redshift and high spatial resolution (~ 
0.01") observations with ALMA will allow us to 
observe redshifted high-J CO lines, J > 9 (see 
ii5.1l for details). For local sources, these lines are 
difficult to detect from the ground. They provide 
important diagnostics for th e power and impact 
of th e accreting black hole (JSpaans & Mciieriny 
20081 ). We use the X-ray dominated region (XDR) 
code of Meiierink fc Spaand ( 2005f) . an improved 
version of the above mentioned XDR-model of 
Malonev et al. ( 19961 ). and solve the ID-radiation 
transport equation for the CO fiuxes through the 
inferred cloud complex. This code has been used 
to cal culate detailed diagnos tics for XDR observa- 
tions ( Meiierink et al.|[2007l ) and was successfully 
applied to infrared and (sub-)milli meter observa- 



tions of luminous infrared galaxies (jLoenen et al 
20081 ). 
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Fig. 5. — The X-ray chemistry in a system with 
X-ray flux of 1 erg s~^ cm~^ impinging on the 
cloud. The adopted density is 10^ cm"'^. Top: 
The abundances of different species as a function 
of column density. Middle: The low- J CO lines as 
a function of column density. Bottom: The high- 
J CO lines as a function of column density. 
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We first calculate the chemical abundances and 
the expe cted CO line i n tensit ies for the model 
used by lOalliano et al. ( 2003 ) as discussed in 
§ 14.11 While they adopted a column density of 
10^^ cm~^, we are extending the calculation to 
larger columns of 10^'* cm~^, corresponding to the 
typical column needed to absorb 5 keV photons. 
This allows us to generalize these results to sys- 
tems with larger column densities and Compton 
thick material. The results of the calculations are 
shown in Fig. |4l We note that the fiducial gas 
density of 10^ cm~^ has little impact on our re- 
1 suits, unless it drops to below lO''-^ cm~^. Above 
this limit, the XDR properties are determined 
^ by the ratio of X-ray flux to gas density. In 
addition, we have checked that even if the en- 
tire disk of ~ 100 pc was at low densities of 
10^ cm~^, the high- J CO line intensities would 



still be 



10 erg s cm 



sr 



as the low- 



density gas extends over much larger spatial scales 
and thus reaches larger or comparable column den- 
sities. 

The strong X-ray flux of ~ 170 erg s^^ cm"^ 
in NGC 1068 suffices to make the gas essen- 
tially atomic and leads to high temperatures of 
~ 3000 K, as well as relatively low CO abun- 
dances of the order 10~^. However, the CO in- 
tensity is high, due to the strong excitation in the 
hot gas. For a column of 10^^ cm~^, our results 
ap pear of the sanie ma gnitude as in the model 
of Galliano et al. 1 (|2003h . For larger columns, the 
temperature gradually decreases, the gas becomes 
molecular and CO gets more abundant, and we 
find intensities of the order 10~^ erg s^^ cni^^ sr~^ 
in the high- J CO lines. 

We need to account for the fact that many sys- 
tems may accrete considerably below the Edding- 
ton limit as discussed above. Also, we want to 
consider quasars with somewhat smaller black hole 
masses of the order 10^ Mq, as their number den- 
sity is larger by an order of magnitude (see Fig.[T]). 
Therefore, we calculate models with less X-ray flux 
as well. An extreme case with ^ 1 erg s~^ cm~^ 
is shown in Fig. [5l In this model, we find lower 
temperatures of ~ 70 K, a large fraction of molec- 
ular gas and CO abundances of the order 10^^. 
While the lower temperature tends to decrease the 
CO line intensities, they are still enhanced due 
to the larger CO abundance. Above a column of 
10^^ cm~^, the intensities increase rather slowly 
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Fig. 6. — The X-ray chemistry in a system with 
X-ray flux of 10 erg cm~^ impinging on the 
cloud. The adopted density is 10^ cm^'^. Top: 
The abundances of different species as a function 
of column density. Middle: The low- J CO lines as 
a function of column density. Bottom: The high- 
J CO lines as a function of column density. 



as the lines become optically thick. 

As an intermediate scenario, we consider an 
X-ray flux of ^ 10 erg s~^ cm~^ (see Fig. [S]). 
In this model, the temperature is increased to 
~ 100 K. The CO abundance is initially of the 
order 3 x 10~^ and increases to ~ lO"'' for larger 
columns. For columns less than 10^^ cm^^, the 
intensities are thus reduced by about an order of 
magnitude compared to the previous case, while 
they are increased by an order of magnitude for 
larger columns. 

From an observational perspective, it is highly 
interesting that even X-ray fluxes of ~ 1 erg s~^ cm 
lead to significant CO line intensities. At a 
column of 10^^ cm^^, the high-J CO line in- 
tensity in all models was of the order lO^"' — 
10"'^ erg s^^ cm^^ sr^^, and it may increase up 
to values of lO^'^ — 10^^ erg s^^ cm~^ sr~^ for 
columns of the order 10^** cni~^. 

5. High-redshift observations v^rith ALMA 
and JWST 



freq. [GHz] 


^rcs [" ] 


Si [mJy] 


Sc [mJy] 


^bcam 


84- 116 


0.034 


8.9 


0.06 


56 


125 - 169 


0.023 


9.1 


0.07 


48 


163 - 211 


0.018 


150 


1.3 


35 


211 - 275 


0.014 


13 


0.14 


27 


275 - 373 


0.011 


21 


0.25 


18 


385 - 500 


0.008 


63 


0.86 


12 


602 - 720 


0.005 


80 


1.3 
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Table 1: Frequency range, angular resolution 6'rcs, 
line sensitivity Si and continuum sensitivity Sc 
in mJy for one hour of integration time and pri- 
mary beam size 6'bcam for bands 3 — 9 of ALMA. 3 
more bands might be added in the future, band 1 
around 40 GHz, band 2 around 80 GHz and band 
10 around 920 GHz, which will have similar prop- 
erties as the neighbouring bands. 

Future observations with ALMA and JWST 
provide a new opportunity to find the first quasars 
at high redshift. Their unprecedented sensitiv- 
ity permits observations of less luminous sources 
and at higher redshift and we can therefore find 
quasars with black hole masses of 10^ — 10^ Mq, 
which may be the progenitors of the known sources 
with ^ 10^ Mq black holes. In addition, their 
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high spatial and spectral resolution can provide 


Observable 


A [/im] 


if [mjy] 


Redshift 


detailed information about the host galaxies of 


high- J CO lines 


~ 300 


3-30 


z > 5 


those bright high-redshift quasars that we already 


[0 I] 3p^ ^3 


63.2 


- 20 


z > 5.7 


know. In this section, we predict the main observ- 


[0 III] ^3 


51.8 


- 20 


z > 3.8 


ables for ALMA and JWST, provide estimates for 


[N II] ^'■^ Pi 


121.9 


- 7.6 


z > 2.5 


the number of high-redshift sources and discuss 


[0 1] 3po^3p^ 


145.5 


- 3.5 


z > 1.9 


also the importance of follow-up observations on 


[CII] ^P:i,2^^Pi/2 


157.7 


- 69 


z > 1.7 


known high-redshift quasars. 


[S III] 3p^ _^3 p^ 


33.5 


- 1.5 


z > 11.5 




[Sill] ^P^,2^^Pl,2 


34.8 


- 2.6 


z > 11.1 


5.1. Observables for ALMA 


dust continuum 


- 500 


- 0.05 


z > 



In the frequency range of ALMA, it is possi- 
ble to observe redshifted emission from high- J CO 
lines and many fine-structure lines. The main 
properties of ALMA are summarized in Table [1] 
For one hour of integration time, ALMA can de- 
tect line fluxes of ~ 10 mJy at 3a. This is partic- 
ularly interesting with respect to the high-J CO 
intensities calculated in § 14.21 At redshift 2 = 8, 
a resolution of ^ 0.01" corresponds to a physical 
scale of ^ 49 pc. A source of such size would fill a 
solid angle of ~ 7 x 10~^^ sr. With a fiducial ve- 
locity dispersion of 20 km/s and a line intensity of 
10"'^ erg s~^ cm~2 sr~^, this corresponds to a flux 
of ^ 30 mJy, which is therefore detectable in less 
than one hour. In the calculations above, we found 
lower intensities only for models with very high X- 
ray fluxes and small columns (see Fig. 's [4] to [6]). 
Even in this case, we still expect a significant flux 
of '--^ 3 mJy. 

As discussed by iFinkelstein et aL (2008), there 
are further contributions from dust emission which 
may be relevant for ALMA as well, because the 
continuum sensitivity is hi gher by two orders o f 
magnitude. For NGC 1068. lBrauher et al.l (|2008h 
report a flux of 85.04 Jy at 25 /im and 224 Jy at 
100 fim. For z > 5, such fluxes would be red- 
shifted into the ALMA bands. An interpolation 
to 500 /Ltm yields a flux of ~ 700 Jy, which would 
be redshifted into ALMA band 3 for z ~ 5. For 
the same emission at this redshift, we would there- 
fore expect a flux of 0.06 mJy, which matches with 
the band sensitivity. For objects at z = 8, we 
need to consider emission at ^ 330 /xm which is 
redshifted into band 3, where the extrapolation 
yields ~ 520 Jy. This yields an expected flux of 
~ 0.05 mJy, which is detectable within a few hours 
of integration time. Such emission may be even 
stronger if the dust is heated further by ongoing 
starbursts. 

In addition to the high-J CO lines, it seems 



Table 2: The main observables for ALMA, as dis- 
cussed in § 15.11 The estimated flux (p for the high- 
J CO lines assumes that one typical resolution el- 
ement of ~ 0.01" is comparable to the source size, 
while the other fluxes correspond to the fluxes ob- 
served for NGC 1068 if it were placed at z 8 
(for lines only observable at z > 10, we placed it 
at z = 12). The rest frame wavelength is given as 
A. The line and continuum sensitivities are given 
in Table □ 

likely that other lines can be detected with 
ALMA as well For o ur fiducial AGN NGC 1068, 
ISpinoglio et al. ( 2005 ) measured a number of dif- 
ferent fine-structure lines with the Long Wave- 
length Spectrometer (LWS) on board the Infrared 
Space Observatory (ISO). From their sample, we 
have identified several lines that are interesting for 
observations at high redshift. Assuming the same 
source at larger redshift, wc find typical fiuxes 
that could be detected with ALMA within one 
to several hours (see Table [5]) . The high angular 
resolution of ALMA may even allow to resolve 
and detect the higher fiuxes from the inner core. 
The detection of the O and C"*" fine-structure lines 
will be particularly interesting, as the sum of the 
[O I] 63 /im, the [C II] 158 /xm and the CO lines 
is close to the total coohng rate of the gas, given 
that above 10~^ solar metallicity the cooling via 
H2 rotational lines is only a small contribution. 

Since at high densities in XDRs the thermal 
balance establishes rapidly, the total cooling fiux 
equals the heating fiux. This allows one to in- 
fer the impinging X-ray flux, while the sum of 
the fine-structure and CO lines divided by the 
heating efficiency is always a measure of the X- 
ray fiux that has been transformed into heat. 
The heating efficiency of ~ 10% is largely in- 
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dependen t of energy and the spectral shape of 



the flux ( Dalgarno et al. 19991 ). For a spatially 



resolved source, the linear size of the cooling 
gas can be obtained from the spectroscopic red- 
shift. The corresponding X-ray flux on that size 
scale can then be converted into a luminosity. 
The observed profiles of the line width may be 
used to derive further estimates for the black 



hole mass and the accretion rate (|Kaspi et al . 
2OOOI: iMcLure fc Jarvid l2002t IVesterga.ardl 1200? 



Warner et al.ll2003l : iDietrich fc Hamannll2004D . In 
addition, a detection of multiple CO lines provides 
the density and column density distribution as a 
function of distance from the central engine. 

Of course, the metallicity is also a caveat re- 
garding some of the observables discussed here. 
While this work assumes that it is close to the so- 
lar value, the situation for lo w-metallicity gas has 
been explored previously by ISpaans fc Meiierink 
(|2008l ). For metaUicities Z > 10"^ Zq, they find 
intensities in the range 10^^—10^^ erg cm~^ sr" 
For even lower metaUicities, the CO lines are hard 
to detect, but emission from HeH+ and H+ be- 
comes more prominent. Particularly interesting 
at very low metaUicities is also the strong emis- 
sion in H2, which will be redshifted into ALMA 
band 10 for sources at z > 10 if this band is re- 
alized. While such low metaUicities seem unlikely 
at z ~ 8, it may become more relevant at these 
higher redshift. 

The choice of metallicity may be guided by 
other observables like the amount of dust emis- 
sion, the detection of a strong starburst by JWST 
as well as theoretial considerations regarding the 
observability of CO lines in low-metallicity gas 
(jSpaans fc Meiierink|[200a ). Such estimates may 
be refined further employing follow-up observa- 
tions with other telescopes that search for the UV 
N'*+ and C"^+ lines. Their line ratio provides an 
important diagnostic for the metallicity and has 
been used to derive t he metallicity in the k nown 
high-redshift quasars ( Pentericci et al.ll2002f ). 



5.2. Observables for JWST 

In the frequency range of JWST, it is possible to 
detect redshifted stellar light, the AGN continuum 
and dust. While such detections will not allow to 
identify the object unambigously as a quasar, they 
could be used to provide a source list for more de- 
tailed follow-up observations with ALMA, which 
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Fig. 7. — The observed spectrum of a starburst 
with 10 Mq yr~^ with an age of 5 Myr located at 
different redshifts. This is compared to the NIR- 
Cam sensitivities for Su detection after an integra- 
tion time of one day. For z < 10, such starbursts 
can be detected within a few hours. 



would allow us to do this ident ification through the 
shape of the CO line SED (jSpaans fc Meiierink 
20081 ). In addition, they provide complementary 



information about star formation in these systems. 
Quasars with unobscured broad-line regions may 
also be visible in further lines. In this section, 
we will focus on the Near Infrared Camera (NIR- 
Cam), which will provide high sensitivity of the 
order a few nJy for a 3a detection at 0.7 — 4.44 ^m 
wavelengths after an integration time of 10"* s (see 
Table [3D and a field of view of (2.16')^- It will 
consist of 4096^ pixels, with an angular resolution 
of 0.0317" per pixel for the short wavelengths and 
0.0648" per pixel for the long wavelengths. 

As discussed in ii l4.1[ we expect a star formation 
rate of ~ 10 Mq yr"! in an object like NGC 1068. 
To calculate the corresponding spectrum, we em- 
ploy the evolutionary synthesis code Starburst99r| 
IjLeitherer eraIlll999l : IVazauez fc Leithereill2005fK 
We adopt their standard IMF with an exponent of 
1.3 in the range 0.1 — 0.5 Mq, and 2.3 in the range 
0.5 — 100 Mq and take the spectrum after a simu- 
lation time of 5 Myr, which corresponds to the age 
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•'http: / /www. stsci.edu/scicncc/starburst99/ 
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of the circumnuclear starburst ring in this object. 
We then calculate the expected flux for this star- 
burst at different redshifts. The results are shown 
in Fig. [7] and compared to the NlRCam sensitiv- 
ity. For an integration time of one day, we find 
that such a starburst can be detected with Str at 
redshift z = 15. For redshifts z < 10, even a few 
hours are sufficient for detection. 

For comparison, we also consider a model with 
a star formation rate of ^ 1 Mq yr^^ (see Fig. [5]). 
In this case, an integration time of one day is still 
sufficient to find a starburst at 2 ~ 5. A more ex- 
tended integration time of 10 days even suffices to 
find such starbursts down to z < 10. This is still 
rather short compared to proposed surveys search- 
ing for light from the first galaxi es, which require 



integ ration times of 0.4 years ([Windhorst et al 
20061) . 



For the AGN continuum. ISpinoglio et al.l(|2005h 
present a model that peaks at 0.1 — 0.3 eV with a 
flux of ~ 0.01 Jy. Such radiation will be redshifted 
into the frequency range of JWST for z > 5. We 
would expect a flux of ~ 0.9 nJy at z = 5 and 
^ 0.3 nJy at z = 8. The first would be detectable 
after an integration time of about three days, and 
the latter within a time of about 15 days. In prac- 
tice, this radiation component will of course over- 
lap with the radiation from the starburst and just 
leads to a mild increase in the total radiation bud- 
get. 



Wavelength 


Sensitivity 


0.70 urn 


7.24 nJy 


0.90 Atm 


4.62 nJy 


1.15 ijm 


3.70 nJy 


1.50 iim. 


3.41 nJy 


2.00 fj.m 


2.83 nJy 


2.77 Atm 


3.81 nJy 


3.22 fj,m 


2.75 nJy 


3.56 /^m 


4.08 nJy 


4.44 (Um 


6.18 nJy 



Table 3: The NIRCam sensitivity of JWST for a 
detection at 3cr after an integration time of 10"* s. 
These values assume that the wide filter is used, 
which is appropriate for continuum fluxes as con- 
sidered here. 

Additional contributions arc also present from 
redshifted Lyman a radiation, which are par- 
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Fig. 8. — The observed spectrum of a starburst 
with 1 Mq yr~^ with an age of 5 Myr located at 
different redshifts. This is compared to the NIR- 
Cam sensitivities for 3a detection after an integra- 
tion time of 10 days. This allows 3a detection of 
such starbursts for z < 10. 



ticularly enhanced in the presence of a star- 
burst or a supermassive black hole. For ex- 
ample in a halo con taining Mh — 10^^ Mq, 
Dijkstra et al. (j2006al |bl) find Lyman a luminosi- 
ties of ~ lO'^^ erg s~^. Considering such a source 
at high redshift and spreading this energy over a 
frequency range of ^ 5 x 10^ GHz, correspond- 
ing to the width of the detector, yields a flux of 
^ 8 nJy at z = 5, and ~ 2 nJy at 2 = 8. As they 
find the luminosity to scale with MfJ^ , this contri- 
bution becomes less significant in smaller systems. 
Such a flux could be detected within several hours. 

If in addition the Mid Infrared Instrument 
(MIRI) camera is employed, which has a field 
of 1.88' X 1.27' and observes at 5 — 28/im with 
sensitivities ranging from 60 — 8000 nJy in 10'* s 
for a 3(7 detection and with R ~ 3000 an intrin- 
sic line width of ~ 100 km s^*, one may aim at 



detecting the H a line at 656.3 nm. This would 
allow us to derive black hole masses and accretion 
rates from the line width, and is complementary 
to to a derivation of the X-ray luminosity with 
ALMA (iKaspi et al.ll2000HMcLure fc Jarvidl2002l: 
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Vestergaardl2002l: I Warner et al.ll2003l : lDietrich fc HamannI 
20041 ). iKaspi et all (|2000l ) give hght curves of the 
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H a and H /3 fluxes for a sample of quasars. We 
have estimated the fluxes obtained from the first 
five sources in their sample, assuming these ob- 
jects would be shifted to redshift z ~ 8. The 
H a line would then be redshifted to ~ 5.9 //m, 
where the MIRI camera has a sensitivity of 
60 nJy for 3ct detection in 10^ s. As we find 
H a fluxes between 100 — 700 nJy, a detection 
is possible within a few hours or less. The ve- 
locity dispersion could be measured with the 
MIRI spectrograph, as the corresponding fluxes 
of 2 X 10""'^^ — 10~^^ erg s~^ cm~^ are well above 
the sensitivity of 2 x 10^^^ erg s~^ cm"-^ for 3cr 
detection in 10^ s. 

For the H /3 line at 486.1 nm, we find that 
it would be redshifted to ~ 4.4 /im, which is 
still in the range of NlRCam, with a sensitivity 
of 6.18 nJy. Our estimates yield fluxes between 
2 — 20 nJy, detectable within some hours up to a 
day. At higher redshift, this will also fall into the 
MIRI frequency range, where it may be detected 
even on a shorter timescale due to the smaller in- 
trinsic linewidth of the instrument. 



Observable 



ip {nJy\ 



estimated from 



l ine, we use the z = 3.4 sample of iDietrich et al 
(|2003bl ). We flnd similarly high fluxes of the order 
20 — 100 /iJy. They may be reduced for less mas- 
sive black holes, but are well accessible for JWST. 

The velocity dispersion of the H /3, Mg+ and 
C'^^ lines can be measured with NIRSpec, that 
provides a sensitivity of 1.5xl0~^^ ergs~^ cm~^, 
while the expected flux is (1 — 10) x 10~^^ erg s~^ cm 
for H /3, (5 - 20) x lO'^^ erg s'^ cm^^ for 0^+ 
and (1 - 8) X 10~" erg s'^ cm^^ ^j. j^g+_ 

In summary, we find that the component from 
the starburst will likely be the main observable for 
JWST, as well as the H a, H /3, Mg+ and 0^+ fines 
in case of unobscured broad-line regions. Even 
the other components alone would be detectable. 
At redshifts z 5, typical integration times may 
reach from a few hours to a day, depending on 
the luminosity of the source. Near z ~ 10, typical 
integration times may vary between one and ten 
days. The main observables are summarized in 
Table H 

5.3. The number of high-redshift sources 
for diflferent observational strategies 



stellar light 
AGN continuum 
Lyman a 
H a 
H /? 

Mg II 279.8 nm 
C IV 154.9 nm 



0.3-3 
0.3 

- 2 

100 - 700 
2-20 

lO-"^ - 3 X 10^ 
2 X 10-* - 10^ 



Davies et al. 2007 
Spinogho et al. 2005 
Dijkstra et al. 2006 
Kaspi et al. 2000 
Kaspi et al. 2000 
Dietrich et al. 1999 
Dietrich et al. 2003 



~ With the considerations above, we can attempt 
to estimate the number of observable sources at 
differe nt redshifts. We adopt a s imilar formal- 
ism as IChoudhurv fc Ferraral (|2007t ) and count the 
number of black holes A^(z, Mbh) with masses 
larger than Mbh in a redshift interval [z, z + Az]. 
This is given as 



Table 4: The main observables for JWST, as dis- 
cussed in § 15.21 We estimate the fiux at z ~ 8 
as described in the text. The NIRCam continuum 
sensitivities are given in Table [31 

In addition, one may search for the Mg+ line 
with a rest frame wavelength of 279.8 nm and 
the line with 154.9 nm, which would be red- 
shifted into the NIRCam frequency range and pro- 
vide complementary information about the veloc- 
ity dispersion near the black hole. We find that, 
if the objects from the sample of IPietrich et al. 

would be shifted to z ~ 8, the expected 
fluxes from the C^^ line in a typical NIRCam 
frequency bin (wide filter) would be of the order 
1 — 3 /xjy, and thus detectable within much less 
than an hour. To estimate the fiuxes in the Mg+ 



iV(z, Af] 



BH) 



z-(-Az ^ /-A/max 

dz' . . .„ / dXogMu 



dz'dfl 



X fa 



dn 



dlogMh 



[Mbh], 



(9) 



where dV dz'~^ dfl^^ denotes the comoving vol- 
ume element per unit redshift per unit solid angle, 
which is given as ({Peebles 1993) 



dV 

dz'dVL 



D\c-. 

dz 



(10) 



In this expression. Da is the angular diameter 
distance, c is the speed of light and dt/dz = 
1/ (iJ(z)(l -I- z)), where H{z) is the expansion rate 
as a function of redshift. The term dn/ dlogMh 
describes the number density of dark matter ha- 
los per unit mass, and the integration range in 
Eq. [S] is chosen such that the minimal halo mass 
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Fig. 9. — The expected number of sources for the 
first observational strategy focusing on high- J CO 
hne emission. We adopt a field of view of (l')^j an 
active quasar fraction of 20% and a redshifter in- 
terval of Az = 0.5 as discussed in the text. Down 
to z ^ 10, a typical ALMA field of view therefore 
has a high chance to find a number of quasars with 
black holes in the mass range of 10^ — 10^ M©. 
Near z = 6, one may even find more massive 
quasars by extending the search to larger volumes. 
Note that the numbers above scale linearly with 
the adopted fraction of active quasars. 



corresponds to a black hole mass of A/bh for our 
model in § [21 The upper halo mass is chosen cor- 
responding to a black hole mass of 10^ Mq. The 
parameter fa is the fraction of active quasars, as 
introduced in § [3l 

As our reference system, we have adopted the 
known AGN NGC 1068, which has a black hole 
of ^ 10^ Mq and a starburst component with 
~ 10 Mq yr~^. As discussed above, even less ac- 
tive quasars provide a sufficient amount of CO flux 
for detection, and the same holds for smaller star 
formation rates. We will therefore assume that 
typical galaxies containing a 10^ Af© black hole 
should be observable. To be more conservative, 
we will also show the corresponding numbers for 
more massive black holes. 

In the following, we will distinguish between 
two different observational strategies. The first 
aims at finding high- J CO lines with ALMA, 
which will allow us to una mbigously identify a 
given source as a quasar (jSpaans fc Meiierink 
20081 ). The second aims at finding high-redshift 
sources in general and focuses mainly on contin- 



Fig. 10. — The expected number of sources for the 
second observational strategy focusing on contin- 
uum radiation mostly from stellar light. We adopt 
a field of view of (2.16')^ and an active quasar frac- 
tion of 100%, as the relevant observables do not 
require quasar activity. We explore a redshift in- 
terval of Az = 2. In fact, this yields a high number 
of sources that are potentially accessible to JWST 
within the NIRCam field of view. For redshifts 
near 2 = 5, typical detection timescales are one 
day or less, while observations near z = 10 require 
integration times between one and ten days. This 
is within reach of deep surveys. 



uum emission from dust or stellar light. In fact, 
once such sources are found, they can be examined 
in more detail by other observatories like the Keck 
telescope or new generations of X-ray telescopes. 

For the first observational strategy focusing on 
line observations, ALMA provides a band with of 
8 GHz. In general, one may scan a redshift inter- 
val of Az/(1 + z) — Af/i/, where Ai/ refers to the 
difference with respect to the observed frequency 
v. For 90 GHz and z ~ 5, this corresponds to 
Az — 0.5. This becomes even larger when going 
to higher redshifts, and smaller for observations 
at higher frequencies, such that we slightly under- 
estimate the sources at high redshift. We adopt 
this value here as a typical mean. Based on the 
results in § [3l we assume an active quasar fraction 
of /a ^ 20%. Our results can however be directly 
scaled to other values. We consider a field of view 
of (1')^: even though ALMA can probably scan 
an angular region of (3')^ within a few days. The 
result is given in Fig. [9l 

We find that down to z 10 there are good 
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chances to detect several sources with black hole 
masses in the range 10*^ - 10^ Mq within the 
ALMA beam. For an extended survey that probes 
an angular field of ~ (3')^j one may even find 
quasars with Mbh > 10^ Mq, and at the same 
time find enough sources with smaller masses to 
allow a statistical analysis. At redshifts beyond 
10, we find an exponential decrease in the number 
of high-redshift quasars, though we expect that 
ALMA could still find some of them. Searching 
for such systems should be a high priority for fu- 
ture observations, as they are a direct probe for 
the origin of the first superniassive black holes. 

For the second strategy focusing on continuum 
radiation, we are not restricted to active black 
holes and work with fa = 1- This assumes that 
such galaxies have a star formation rate of at least 
1 Mq yr~^. However, our results can simply be 
rescaled for different star formation rates. For con- 
tinuum observations, the bandwidth is no longer a 
restriction for the redshift range, and we consider a 
generic interval of Az ~ 2. We adopt the NIRCam 
field of view of (2.16')^. The result is then given in 
Fig. [ini While we cannot be sure that the detected 
objects are active quasars, it is always possible to 
address this question with follow-up observations 
involving other telescopes. In the NIRCam field 
of view, we expect 10^ — 10'^ observable sources at 
redshifts 5 — 10, and 10 — 100 at redshifts beyond. 
In addition, we expect a few very bright sources 
that could be associated with 10^ Mq black holes 
within the field of view. They can be found with 
integration times of a few hours, but beyond that 
a deep field is required to see smaller galaxies and 
probe the early epoch of galaxy formation in de- 
tail. 

5.4. What we will learn about the bright- 
est quasars 

In addition to finding less massive quasars 
at higher redshift, instruments like ALMA and 
JWST can be used to get more detailed informa- 
tion about those sources that we already know. 
At redshift z — 6, a typical angular resolution 
of ~ 0.01' corresponds to a spatial resolution of 
~ 60 pc. It is thus possible to resolve the quasar 
host galaxy in high- J CO emission and dust with 
ALMA. As these sources are very bright in X- 
ray emission, the required fiux to density ratio 
that produces significant CO emission is shifted 



to larger scales further away from the central black 
hole and fills a larger region in space. With de- 
tailed diagnostics for the CO lines, it is thus pos- 
sible to probe temperature and density in these 
systems with spatial resolution. With a spectral 
resolution of up to 0.01 km s~^, the dynamics of 
the gas can be probed in detail, in particular with 
respect to typical infall velocities as well as winds 
and jets. 

In addition, JWST will provide complementary 
information about stellar light, the AGN contin- 
uum and Lyman a emission. The latter could be 
clearly discriminated from the continuum contri- 
bution using specific narrow filters tuned to the 
wavelength corresponding to the redshifted Lyman 
a line. 

6. Discussion and outlook 

In this paper, we have presented a model for the 
quasar population at high redshift, based on the 
known correlations between the black hole mass 
and its host galaxy. We have derived constraints 
on the fraction of active quasars from the SXRB 
and discussed the expected observables. We have 
based our discussion on using NGC 1068 as a tem- 
plate AGN, and calculated its observability at high 
redshift. Our focus lies in particular on high- J CO 
lines, fine-structure lines of C"*" and O, dust emis- 
sion, stellar light, the AGN continuum and Lyman 
a radiation. Above a metallicity of 10~^ solar, 
CO, 0+ and O are the main coolants, and their 
detection can be used to derive an estimate for 
the heating rate and the impinging X-ray flux. A 
profile of the observed line width further allows 
to derive estimates for the black hole mass and 
probes the mass - luminosity relation at high red- 
shift. We expect that galaxies containing a black 
hole in the mass range 10^ — 10^ Mq can be ob- 
served at high redshift with ALMA and JWST. 
We predict 1 — 10 high-redshift sources that can 
be found via molecular lines in an ALMA field of 
view of (1')^, and 100 — 1000 sources to be found 
by JWST using NIRCam and focusing on stellar 
emission. 

One uncertainty in this discussion comes from 
the question whether the Maggorian relation is es- 
tablished already at high redshift. This seems fair 
in light of supersolar metallicities in high-redshift 
quasars. On the other hand, if our model over- 
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predicts the number of high-redshift black holes, 
the X-ray constraints would be alleviated further, 
and a larger fraction of the black hole population 
could be active. In this sense, the number of active 
quasars might still be of the same order. Indeed, 
unless the population of black holes between red- 
shifts 5 and 6 has accreted mass with high super- 
Eddington rates, there must be progenitors with 
masses of - 10^ - 10^ Mq at higher redshift. It 
is a secondary issue for our purposes whether the 
black hole would be directly visible in X-rays, or 
if the system is a hidden AGN and only visible in 
the excited high- J CO lines. 

In about three years, ALMA will become opera- 
tional with initially 16 antennas. This will provide 
sufficient sensitivity to detect high-J CO lines, 
fine-structure lines and dust continuum emis- 
sion from the brightest high-redshift quasars that 
are already known. As these systems are large, 
the high angular resolution of ALMA will pro- 
vide detailed spatial information about the inner 
structure of such systems, and detailed diagnos- 
tics for the CO lines allow to infer the density 
and temperature dist ribution in the host galaxy 



( Meiierink et al.ll2007n 



Once ALMA has the full set of 50 antennas, it 
reaches even higher sensitivity and it becomes pos- 
sible to search for AGNs with less massive black 
holes. For such deep field observations, we rec- 
ommend to focus initially on modest redshifts of 
z ~ 5 — 6, where we expect higher fluxes, more 
sources, and where one may be more confident 
that the Magorrian relation is well-established. 
This provides an important test to confirm the ob- 
servability of high-redshift quasars, as well as im- 
portant information about the progenitors of the 
most massive quasars in the universe. Such studies 
can then gradually be extended to higher redshifts 
and will provide a direct test for structure forma- 
tion models in the high-redshift universe and of 
the Magorrian relation at different redshifts. 

JWST can aid such surveys by providing ad- 
ditional sources that may be detected in stellar 
light, the AGN continuum or Lyman a emission. 
In this case, ALMA does not need to glance into 
the dark sky, but can do pointed follow-ups of the 
known sources. This will speed up high-redshift 
surveys, and provide complementary information 
at different frequencies. It is particularly impor- 
tant at very high redshift beyond z 10, where 



the probability of finding sources in the ALMA 
field of view decreases rapidly. The combination 
of these instruments can therefore truely detect 
the first quasars in the universe. 
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